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>■ ' Abstract 

' In cascade inflation and some other string inflation models, collisions of mobile 



branes with other branes or orbifold planes occur and lead to interesting cosmologi- 
cal signatures. The fundamental M/string-theory description of these collisions is still 
lacking but it is clear that the inflaton looses part of its energy to some form of brane 
. matter, e.g. a component of tensionless strings. In the absence of a fundamental de- 

■ scription, we assume a general barotropic fluid on the brane, which absorbs part of the 

^ ' inflaton's energy. The fluid is modeled by a scalar with a suitable exponential potential 

' to arrive at a full-fledged field theory model. We study numerically the impact of the 

energy transfer from the inflaton to the scalar on curvature and isocurvature pertur- 
bations and demonstrate explicitly that the curvature power spectrum gets modulated 
by oscillations which damp away toward smaller scales. Even though, the contribution 
of isocurvature perturbations decays toward the end of inflation, they induce curvature 
perturbations on scales that exit the horizon before the collision. We consider cases 
where the scalar behaves like radiation, matter or a web of cosmic strings and discuss 
the differences in the resulting power spectra. 
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1 Introduction 



There has been a lot of activity over the past few years to derive inflation from string 
theory. Most prominent among the many approaches has been brane inflation, in which 
the coordinates of some mobile branes on the internal compactification manifold provide 
the inflatons in the effective four- dimensional low energy theory [Ij. While initial efforts 
focussed, for simplicity, on compactifications with a single mobile brane (for recent updates 
on single brane inflations, see e.g. [2]- [10]) general compactifications with fluxes naturally 
possess several, often many, mobile branes to satisfy tadpole cancellation equations p!T]-[T6]. 
Multi brane inflation models lead to multi field inflation cosmologies in the low energy theory, 
which as a result have been actively researched recently (see e.g. [IZ]-[3S])- The multitude 
of inflaton fields can be actually seen as a physical blessing rather than a technical curse. 

The multi brane inflation models split into two classes: those with brane-brane interac- 
tions which steer the multi brane system towards a stable dynamical attractor and those 
with brane-brane interactions for which no such attractor exists. The first class allows to 
replace the multi field cosmologies by cosmologies involving fewer fields, in the extreme case 
just a single inflaton field, once the system evolves along the attractor. An example for the 
first class is a multi brane system with exponential interactions among the branes, caused 
by non-perturbative interactions between classically non-interacting branes pULj- As a result 
of the stable attractor evolution this class more readily allows for an analytical treatment 
and moreover offers a parametric way to achieve inflation. Namely by choosing the number 
of participating branes large enough, one achieves a parametric suppression of the slow-roll 
parameters - a phenomenon which is known as assisted inflation [10] . 

The second class, on the other hand, requires an enormous fine-tuning for the position of 
each mobile brane in order to yield an ordered, e.g. equidistant, brane configuration, which 
admits an analytical treatment. An example of this class would be a system consisting 
of many brane-antibrane pairs with mutual Coulomb interactions. Full-fledged numerical 
studies for the cosmic evolution of this second class still need to be performed. Nevertheless, 
one might expect that this second class won't lead to inflation. The reason is that even if 
one starts with a multi brane configuration suitable for inflation (which, as said, requires 
a tremendous fine-tuning for all brane positions as the system lacks an attractor), small 
perturbations won't die away. They will lead the system to a complicated multi brane 
dynamics which moves the branes around on the internal space in an unconcerted way, 
precluding the realization of assisted inflation. The upshot might be that among the many 
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possible multi brane models, the models of the first class are the ones capable of generating 
inflation, whereas the models of the second class have to be dismissed. 

In this paper we focus on a new phenomenon which arises in multi brane inflation models: 
the possibility of cascade inflation phases [13]. These arise when some of the mobile branes 
collide successively with fixed branes or orbifold fixed planes and thus no longer participate 
in the inflation process. The parameters of the effective four-dimensional inflationary po- 
tential vary in these collisions and as a result the inflationary potential acquires features. 
A characteristic consequence of such features is a damped oscillatory behavior of the power 
spectrum of density perturbations. This could lead to important observational clues about 
the underlying time variation of the inflationary potential. 

To date, we have no complete fundamental M/string-theory description of such collisions. 
Therefore, we have to model the collision in field theory as an energy transfer (ET) from 
the inflaton to the brane. For instance, the brane could be modeled by a perfect fluid 
with a coupling to the inflaton which allows to transfer the inflaton's energy in discrete 
steps. More adapt to a field theory description based on a Lagrangean, and the route 
we follow here, is to model the perfect fluid by a scalar field with appropriate potential. 
Our goal in this paper is to study the potentially observable imprints of this ET on the 
resulting curvature and isocurvature perturbation spectrum. Chapter [2] describes the M- 
theory motivation and background for our analysis, reviewing cascade inflation. Chapter [3] 
introduces our two scalar field theory, which models the ET in brane collisions such as arising 
in cascade inflation. Chapter H] presents the curvature and isocurvature perturbations for 
this model. We find that the contribution of isocurvature perturbations decays toward the 
end of inflation. Nonetheless, they induce curvature perturbations on scales that exit the 
horizon before the collision. The curvature power spectrum gets modulated by oscillations 
which damp away toward smaller scales. Finally, technical details on how to calculate the 
perturbations in our two-field model are presented in the Appendix. 

2 Cascade Inflation 

The goal of this paper is to study perturbations generated in a two scalar field inflation model 
with exponential potentials. Such potentials arise in fundamental physics from M-theory 
multi brane inflation and cascade inflation [TB] through non-perturbative instanton in- 
teractions. Since the two scalar field inflation model might also be used to provide some first 
insight into the cosmological impact of the ET arising when a bulk M5-brane collides with 
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the boundary in cascade inflation, we will now briefly recall cascade inflation as one of the 
motivations for our study. 

Cascade inflation arises from heterotic M-theory [17], |18] compactified on X x S^/Z2, 
where X is a Calabi-Yau threefold. Tadpole cancellation requires generically the presences 
of M5-branes in the background. These fill the four-dimensional spacetime and wrap a 
holomorphic two-cycle on X of genus zero. As BPS objects, the M5-branes interact only 
non-perturbatively via open M2-instantons which stretch along S^/Z2 and wrap the same 
two-cycle^]. The effective four- dimensional A/" = 1 supergravity resulting from this compact- 
ification contains as scalar components of complex chiral superfields the M5-brane position 
moduli Yi, i = 1, . . . , N, the Calabi-Yau X volume modulus S and the S^/Z2 orbifold-size 
modulus T. Twice their real parts are denoted by yi, s and t. Moreover, it's convenient to 
define i/ = (Eti 2/')'^'- One can show that the M-theory dynamics leads to an equidistant 
distribution of the M5-branes along S^/Z2 [H]. For such a distribution 

Re(y,+i-F,)= (^^^Ax (1) 

is independent of the M5-brane counting label i. Consequently, the A^ — 1 inflatons Re(yj+i — 
Yi) can be identified. The multi inflaton system reduces effectively to one with a single 
inflaton Ax. The S^/Z2 interval size is denoted by L. 

In the large volume limit, specified by the inequality st ^ y'^, where supergravity provides 
a reliable description of the dynamics [TT], the potential for the canonically normalized 
inflaton Lp ~ N^^'^Ax becomes 



VN{^) = VNe Vp^mp , (2) 

where 

and d is the Calabi-Yau intersection number. For the detailed derivation from M-theory we 
refer the reader to [HI [13] . The cosmological FRW evolution in this background is given by 
power law inflation [19] with FRW scale factor 



a(t) = aot^^ . (4) 

Inflation sets in when > 1, which can always be achieved when sufficiently many M5- 
branes are present since Pat ~ A^^. Besides this bound on A^ from below, there is also a 



^For simplicity we are assuming a Calabi-Yau threefold with Hodge number /i(^'^^(X) = 1. 



4 



bound on from above which follows from the requirement to work in the large volume 
regime where st ^ y'^ holds and the fact that y grows with A^. For typical parameter 
values one thus finds 20 < < 200 as a constraint on the number of M5-branes [H] . Such 
numbers can easily be accounted for in heterotic M-theory flux compactifications [50l [STl [52] 
where tadpole cancellation equations balance the amount of M5-branes with quantized flux 
numbers. 

The repulsive M2-instanton interactions between the M5-branes cause them to spread 
over the S^/Z2 interval until the two outermost M5-branes hit the boundaries and dissolve 
into them via small instanton transitions [53], [SI], [5S]- This process changes the topological 
data on the boundaries while the number of M5-branes participating in the inflationary bulk 
dynamics drops from to N — 2. The remaining N — 2 bulk M5-branes will continue 
to spread until again the most outermost M5-branes hit the boundaries in a second small 
instanton transition and so on. This evolution, in which the number of M5-branes drops 
successively in discrete steps, defines cascade inflation [13j. In [13j, the analysis neglected 
the energy transferred to the boundaries by the instanton transitions and therefore worked, 
after each instanton transition, with a suitably modified power-law evolution 

'^m(t) ~ 0,m^^^"^ , tm— 1 ^ t ^ t^^ , (5) 

having different A^ dependent parameters pn and Vat for each interval. The cascade inflation 
process terminates when the number of M5-branes in the mth phase, Nm = N — 2m, drops 
below a critical value Nk in the final i^th phase. This critical value Nk is determined by 
the exit condition ptv^ = at which infiation stops and which is dynamically reached from 
larger powers p^^ ^ Pnki ^ — ^ ■ Thus we have a finite number m = 1, . . . ,K oi cascade 
infiation bouts. 

Throughout the whole cascade infiation process the infiaton. Ax, representing the dis- 
tance between neighboring bulk M5-branes, grows continuously. Matching the scale factors 
at the transition times, t^, determines the prefactors 

_ .PN^-PN2.PN2-PN^ PN^_-^-PNm /^N 

U-m — U-iL]^ ■ ■ ■ i I 

where iij = tj/tj. The scale factor, but not the Hubble parameter, is continuous at the 
transition times t^, when the ET to the boundaries is neglected. The onset time of infiation, 
to, is determined by inverting the power-law infiation solution for (p{t) in the initial phase 
and noting that Ax(to)/L <^ 1. The result is 



Inverting the solution for (f{t) gives the transition times 

(t„. - t.)Mp = (PElpfLZ^lA + ± iAzlle-i^--^) M . (8) 

^ ^ k=2 ' 

from which the number of e-foldings, generated during cascade inflation, derives 

iVe = In = X^P^. ln(t^,^_i) . (9) 

^ \ 0/ / ■m=l 

The analysis of [13] neglected the backreaction of the energy, which gets lost to the 
boundaries. The reason for this was the still open question how the small instanton tran- 
sitions should be described dynamically at a fundamental level plus the naive expectation 
that, as long, as the number of M5-branes, having been absorbed by the boundaries, stays 
small compared to those remaining in the bulk and driving inflation, this might be a useful 
approximation. For the type I string, possessing an S'O (32) gauge group, the small S'0(32) 
instanton is nothing but a D5-brane in ten dimensions [5^. When compactified on a K3 
manifold from ten to six dimensions, a six-dimensional instanton arises from the D5-brane 
wrapping the K3 which describes an effective string in six dimensions (as the string theory 
situation is best understood in six dimensions, we stick to this case for this brief discussion). 
At a certain point in moduli space the tension of this string vanishes and a singularity occurs. 
At this point the string's sigma model coupling is of order one, hence the sigma model is 
strongly coupled, whereas the string coupling constant could be arbitrarily small [57]. The 
type I - heterotic duality relates this type I phenomena to a singularity involving the fun- 
damental 5*0(32) heterotic string in six dimensions. Furthermore, invoking T-duality, one 
can relate the phenomenon further to the x E'g, heterotic string in six dimensions [56] , for 
which at the singularity one of the two gauge couplings diverges and the string associated 
with the small instanton in that gauge group becomes tensionless, while the other E^ gauge 
coupling remains finite. 

In heterotic M-theory, compactified on K3xS^/Z2 from eleven down to six dimensions, 
the two E^ gauge groups are geometrically separated and localized on the ends of the S^/Z2 
interval. For these we can first of all have the same instanton in either E^ gauge theory as for 
the weakly coupled E^ x E^ heterotic string with the associated tensionless string when the 
E^ instanton shrinks and produces a singularity. However, a novel situation arises from the 
presence of M5-branes, which are generically needed to cancel anomalies. In six dimensions, 
these M5-branes are located as points on the compactification manifold K3xS^/Z2. An 
M2-brane stretches along S^/Z2 from the M5-brane to each boundary. Such an M2-brane 
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produces a string in the non-compact six dimensions, and generates a tensionless string 
once the M5-brane hits a boundary [56], [57]. In heterotic M-theory compactifications on 
X X S^/Z2 down to four dimensions, the small instantons are described by a torsion free 
sheaf, a singular bundle. The singular torsion free sheaf can then be smoothed out to a 
non-singular holomorphic vector bundle by moving in moduli space |54j . 

To date, unfortunately, no clear fundamental M-theory description of these small instan- 
ton transitions is available, which would fully describe its dynamics, including the produced 
tensionless strings. In what follows, switching back to a four dimensional analysis of the 
ensuing cosmology, we will therefore adopt a quantum field theory description which models 
such a transition by coupling the inflaton ip to another field x (which in heterotic M-theory 
would come from the boundary). To allow for an energy transfer from the inflaton ip to the 
boundary field x, we will introduce a suitable couphng between the two fields. 

3 The Two Field Model 
3.1 The Potential 

As mentioned above, repulsive M2-instanton interactions cause the M5-branes to spread over 
the S^/Z2 interval until the two outermost M5-branes hit the interval's boundaries. The 
ensuing non-perturbative small instanton transition transforms the outermost M5-branes 
into small instantons on the boundaries [53]. This process changes the topological data on 
the boundaries while the number of M5-branes, participating in the inflationary bulk 
dynamics, drops to — 2 at each such transition. The resulting features in the infiaton's 
potential lead to interesting new observational phenomena. These have been derived in [13] 
under neglect of the full backreaction of the energy which gets transferred to the boundaries. 
To get a first clue how this backreaction alters the observational signatures of [13], we will 
now introduce a simple two field model which allows us to study the backreaction effects 
numerically. In our subsequent analysis we will restrict ourselves to a single collision, i.e. a 
single step in the inflaton potential. Furthermore, to obtain a numerically treatable model 
in which we can study the cosmological implications of the ET, our model will be based 
on just two degrees of freedom, the inflaton Lp and a scalar field X) representing the energy 
absorbing boundary. 

The field x is used to model in field theory a variety of boundary barotropic fluids by 
endowing x with a suitable exponential potential, as we now explain. Since we do not know 
the equation of state of the interacting tensionless strings which are being produced in the 
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small instanton transitions and into which the energy of the inflaton is fed, one would like 
to adopt a general barotropic perfect fluid with equation of state 

P = wp, (10) 

which will absorb part of the inflaton's energy. The influence of the ET on the cosmological 
perturbations can then be studied for various values of the parameter w. This way one 
could e.g. model radiation, matter or cosmic strings as possible effective components being 
generated on the brane in the collision. Now in a field theory framework one should rather 
describe the perfect fluid by a suitable field. In fact, it had been shown in |19| that a scalar 
field X with an exponential potential 

V{x) ~ e^V^^ (11) 

leads to a power-law evolution of the scale-factor, a(t) ~ t^, even when | < g < 1. This type 
of evolution would also result from a perfect fluid with equation of state parameter 

w = ^-l. (12) 
3q 

Hence we will model the perfect fluid in a fleld theory description by the scalar x together 
with the above exponential potential. The reader should bear in mind that the energy content 
of the fluids are subdominant with respect to the inflaton's energy and thus the transferred 
energy does not dominate the evolution of the background. Our work thus differs from 
others, such as [59j, in that respect. 

Our potential for the two scalar fleld model 

V{^,x) = Vi^) + W{x,v) (13) 
consists therefore of two different components. First, there is the inflaton potential 




- "(^)-H . (14) 

with ips being the value of the inflaton at the step (brane collision). This type of potential 
results from an M5-brane collision in cascade inflation, see eq. ([2]). The amplitude Vn and 
parameter vary with in a collision. We model this step- like variation using an inflaton 
dependent amplitude and parameter (the index i refers to the initial state before hitting the 
step, whereas the index / refers to the flnal state after hitting the step) 



with a tanh dependence, which interpolates between —1 and 1 and smoothes the step, while 
the parameter Aip governs the smoothed step's width. As said earlier, we will, for simplicity, 
focus our attention on a single step (single brane collision) in the inflaton's potential. 

Second, there is the scalar x which has to absorb a certain amount of energy in the course 
of the collision, as can be seen as follows. Before the collision, there exists only the inflaton, 
(p, with an exponential potential with parameter pi and amplitude Ui. During inflation the 
inflaton's value increases until it approaches the step at a value ips- Here the exponent drops 
to pf and the potential's amplitude to Uj. The difference between initial and final inflaton 
potential energy has to be absorbed by the boundary fluid, i.e. is transferred to the field, x- 
The potential for x should be of the form given in eq. (ITTi) for the reasons explained above. 
However, since energy is transferred to Xi fhe amplitude of its potential will change in time 
and thus depend on ip. Similarly, we can expect the nature of the fluid to change during the 
collision as e.g. a tensionless string component is created which hadn't been there before the 
collision. Therefore, we model the potential for x by 



with inflaton dependent amplitude w{ip) and parameter r{ip). These are again expressed 
in terms of tanh functions 



to describe a smoothed out step. In this way we arrive at a coupling between Lp and x- 

Before encountering the step, the scalar x describes a fluid with equation of state parame- 
ter = g|- — 1. The potential of x is almost zero because its amplitude is almost vanishing. 
This changes quickly during the collision, at which the amphtude raises to Ui — Uf (the 
ET coming from the inflaton) and the equation of state parameter soon adjusts itself at 
w = — 1. In what follows, we will mostly focus on the case where the inflaton's energy 
is transferred to an effective radiation component {q = 1/2) but will also consider matter 
(g = 2/3) or a network of cosmic strings component (g = 1). We concentrate on the radiation 
case first and choose for definiteness the following parameter values 




(17) 




(18) 



(19) 



Pi = 40.138 , 



Pf = 36.598 , 




0.068 , 



(20) 
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Figure 1: The graphs show the evolution of the Hubble parameter and the first slow- roll 
parameter as a function of the number of e-foldings. Around N^, ~ 10 the infiaton's potential 
energy {Ui ~Uj) is transferred to the x field. 



The initial value of x is chosen to be zero, Xi = 0? such that the infiaton's potential energy 
Ui — Uf is transferred to the x field. 

3.2 Cosmological Evolution 

Fig. [T] shows the evolution of the Hubble parameter and the first slow-roll parameter of the 
model as a function of the number of e-foldings. Around — 10 the step is encountered 
(brane collision takes place) and the fraction {Ui — Uf)/Ui of the infiaton's potential energy 
is transferred to the x field. Since the x direction in the total potential V{ip, x) is much 
steeper than the inflaton direction, its energy content redshifts within a few e-foldings and 
the background evolves solely under the influence of the final inflaton potential 



after that. As the slope of the potential increases after the energy of x redshifts, the slow-roll 
parameter, e, settles to a slightly larger value, see the right graph in fig. [H 

The trajectory of and x i^i field space and the evolution of their combination 6, which 
had been defined in eq. (jH]), is graphed in fig. [2l There is a sharp turn in the classical 
trajectory of the fields, when the inflaton, ip, transfers part of its energy to x- As we will 



ifs = 1.477Mp 



(21) 




(22) 
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Figure 2: The left graph shows the evolution of the fields tp and x- They move upwards from 
left to right along the black trajectory, whose very first part coincides with the axis. The 
right graph displays dO/dNf. as a function of the number of e-foldings A^^e- A sharp turn in 
field space at the collision is clearly visible. 

see, with this bent in the trajectory, the curvature perturbations are strongly fed by entropy 
perturbations. We will now calculate the curvature and isocurvature spectra for our energy 
exchanging two field model. 

4 Adiabatic and Isocurvature Perturbations 

For the two-field model described above, we integrate the equations of motion for the cur- 
vature and the isocurvature perturbations , eqs. (1551) and fISBl) . whose derivation is outlined 
in Appendix. 

4.1 Evolution of Curvature and Entropy Perturbations. 

Fig. [3] shows the result for the evolution of curvature and entropy fiuctuations for two modes: 
one that exits the horizon before the ET at \og{k / aoHo) = 7.5, and a second mode which 
exits it after the ET at log(/c/ao-ffo) = 10.5. The subscript denotes the values of the scale 
factor and Hubble parameter at the beginning of infiation. Since we have to deal with two 
independent physical degrees of freedom, we perform the integration twice to account for 
the two independent quantum fiuctuations. In the first run, we assume that Q„ is initially 
in the Bunch-Davies vacuum and that 6s initially vanishes. In this manner, we obtain the 
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log(k/(aoHo))=7.5 




log(k/(aoHo))=10.5 




Figure 3: The left and right graphs show the evolution of IQJJ'^™], \Q'^'^\ and |5sP™| 
as functions of Ne for two comoving wave-numbers that exit the horizon before respectively 
after the ET. 

primordial curvature perturbations Qp'"'™ and the induced isocurvature perturbations Js™*^. 
In the second run we interchange the initial conditions for and 6s to obtain Q™*^ and 
^_gpnm_ rjnj_^g total amplitude of the curvature perturbations is given by: 



|gtot|2 _ |gind|2 ^ |gprim|2 



(23) 



and the amplitude of the total isocurvature perturbations is given by an analogical expression. 
In the single-field case, we perform just one integration. 

Deep inside the Hubble radius, the two perturbations evolve independently in the same 
way, up to a slow overall rotation which practically does not change the amplitudes or the 
correlations. After the Hubble radius crossing approaches the value it would have 

obtained in absence of any isocurvature perturbations. The induced perturbations, Q™*^ and 
^^ind practically negligible inside the Hubble radius. The former can be generated on the 
super-Hubble scales, where its EOM in the slow-roll approximation reads [61], 



(24) 



given that the background trajectory in the field space is sufficiently curved. Also \Q]^'^\ 
for such modes is much smaller than Q^^^™ and consequently we do not see considerable 
enhancement at such scales. 

The behavior of the field valued function dO/dN is displayed in the right graph of fig. [2l 
Before the ET, this function is zero and thus no induced curvature perturbations, IQ"'^], are 
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Figure 4: The left graph shows the adiabatic spectrum vs. log{k/aQHo) for the modes 
that exit the horizon around the decay time. It is assumed that the energy of the decay 
products redshifts as radiation, Uf/Ui = 0.932 and Acf) = 10^'^ Mp. The right graph shows 
the adiabatic power spectrum for the single field case with an inflaton potential having a 
step of equal height. The oscillations in the single field case last much longer than in the 
two field case with ET. 

generated. As the ET from ip to x happens, the sharp turn in the classical trajectory creates 
a spike in dO/dN. This will lead to a considerable increment in the curvature perturbations 
on super-Hubble scales, due to interaction with isocurvature perturbations. As the energy 
of the X" field redshifts, dO/dN becomes zero again and the amplitude \Q^^'^\ becomes frozen. 
This results in an overall enhancement of the amplitude of the curvature modes that exit 
the horizon before the ET. For modes that exit the horizon after the ET, the function \Q^^^\ 
undergoes some modulated oscillations before becoming constant at super-horizon scales. 
As we will see, these oscillations will imprint themselves as modulated oscillations on the 
curvature spectrum. 

4.2 Curvature and Entropy Spectra 

The left plot of fig. H] demonstrates the dependence of the adiabatic spectrum on log(A;/ao-ffo)- 
For comparison, we have also shown the power spectrum for a single inflaton model with 
exponential potential possessing a step of same height. This model is obtained by setting 
the potential for the x field, W{x^({>), and therefore also the coupling to ip, to zero. The 
height of the potential step in the single field case is entirely transferred into kinetic energy 
of the inflaton [601. One point that easily gets noticed by comparing these two graphs is 
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Figure 5: The left graph shows the curvature spectra vs. log(/c/aoi/o) ioT AU/Ui = 0.068 
and AU/Ui = 0.03. The right graph shows the adiabatic power spectra vs. log{k/aoHo) for 
A(p = lO^^Mp (black solid hue) and A(p = lO^^Mp (dashed grey line) 



that the modulated oscillations in the power spectrum for the two field case decay much 
faster than in the single field case. For Uf/U = 0.932, the modulated oscillations last for 
four decades in k for the two field case, whereas in the single field case they continue for 
more than eight decades. Also the amplitude of those modes that exit the horizon before the 
decay is increased by 41%. As mentioned above, this is due to the strong interaction of the 
adiabatic and isocurvature perturbations, when part of the infiaton's energy is transferred 
to X- 

In the slow roll approximation, the EOM of the isocurvature perturbations on the super- 
Hubble scales can be approximated by: 

dSs 



dN, 



-VssSs . (25) 



e 



Since the potential is typically much more curved in the direction orthogonal to the trajec- 
tory in the field space than in the direction along the trajectory. Hence, the isocurvature 
perturbations decay exponentially from their corresponding value at the horizon crossing. 
The shape of the isocurvature spectrum is given in fig. [61 Their corresponding value is 
sub-dominant with respect to the curvature perturbations by a factor which varies between 
10^'^ to 10"^ at different scales. One can see a slight amplification followed by modulated 
oscillation at the scales that leave the horizon around the ET. 
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We now consider the effect of changing 

Af//f/ = (f/, - Uf)/U, , (26) 

keeping all other parameters constant. The left graph of fig. [5] shows the adiabatic spectra 
for Af//f/ = 0.068 and Af//f/ = 0.03. Reducing the energy which gets transferred to the x 
field by a factor of ~ 2.25, decreases the amount of increment in the adiabatic perturbations 
at wave numbers smaller than /cet = a-f^|ET by a factor of ~ 4.2. The amplitude of the mod- 
ulated fluctuations at A; > /cet decreases too, even though the frequency of the oscillations 
remains more or less the same. 

We also considered the effect of changing the parameter Ay? which physically corresponds 
to the decay width of the inflaton. The result is displayed in the right graph of fig. O The 
smaller this parameter is, the faster the energy from inflaton, y9, transfers to x- By decreasing 
this parameter, the adiabatic power spectrum is shifted slightly toward larger scales. This is 
intuitively understandable, as in this case the inflaton energy, A?7, is exchanged faster and 
thus the resulting oscillations start at larger scales. 

So far we have assumed that the decay product x has a potential which makes it redshift 
like radiation, q = 1/2. We now relax this constraint and consider how the perturbations 
will evolve if we choose q appropriate for matter, g = 2/3, or a web of cosmic strings, q = I- 
Fig.(j6j) shows how the curvature spectrum changes once x decays like radiation. For a flxed 
amount of energy transferred to x fleld, increasing q, reduces the amount of ampliflcation 
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Figure 7: The left graph shows the tensor spectra vs. log(A;/aoi/o) for ^U/Ui = 0.932. Sohd 
and dashed hnes respectively represent the single field and two field cases. The right plot 
shows the dependence of the slow- roll parameter e on N^,. 

of curvature spectrum at wave-numbers k < /cet- Thus the least amount of amplification at 
such scales occurs for q = 1. The amplitude and frequency of the oscillations are more or 
less independent from parameter q. 

4.3 Tensor Spectrum 

Finally, we are investigating the tensor spectra of the energy exchanging two field infiation 
model. The left plot of fig. [7] shows the profile of gravity waves for the modes that exit the 
horizon during the ET. For comparison we have also plotted the tensor spectrum for the 
single field case, in which the energy of the step in the ip potential is snatched by the kinetic 
energy of itself. In the single field case the spectrum displays some oscillations for the 
modes that exit the horizon around the ET. This could be understood by comparing the 
variations of e for the two cases. Tensor perturbations satisfy the following equation [62] 

P'l + [k' - ^) P. = 0. (27) 

The quotient a"/ a can be written in terms of the slow-roll parameter e as 

— = 2a^H\2 - e) . (28) 
a 

To understand its implication for the tensor spectrum, we exhibit in the right plot of fig. [7] 
the evolution of e with Ne for the two field and single field cases. In both cases, during the 
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non-slow-roll phases, there are sharp spikes in e. However, for the single field case, the spike 
in e is much greater than in the two field case. This large variation in e in the single field 
case leads to modulated wiggles on the amplitude of its tensor spectrum and explains the 
difference between the single and two field cases seen in the left part of fig. [71 

This result for the single field case should be contrasted with the results of [60], where no 
observable signature in the tensor power spectrum had been observed. However, one should 
note that the energy difference in the step was chosen much smaller in [60] and thus e would 
not exhibit such a sharp spike. 

Furthermore, it can be seen from the left plot in fig. [7] that the amplitude of the tensor 
spectrum decays much faster in the single field than in the two field case. This is because the 
liberated energy from the potential's step transforms in the single field case to the kinetic 
energy of the infiaton which redshifts like a~®. On the contrary, in the two field case the 
released energy is absorbed by the x fi^ld which redshifts much slower. This will cause the 
the Hubble parameter in the single field case to diminish much faster. As the amplitude of 
the tensor perturbation is roughly given by if/27r, this explains the smallness of the tensor 
spectra in the single field case compared to the larger tensor spectra of the two field case 
which we see in the left plot of fig. [71 
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Appendix: Curvature and Isocurvature Perturbations 
in Two-field Infiation 

The calculation of cosmological perturbations in the multi-field infiation is an extensively 
studied topic. Nonetheless, we would like to review the basic notation, results and, in 
particular, the equations of motion for the perturbations that we solve numerically. In this 
section, we shall follow closely the presentation of [H]. 
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A two-scalar-field system coupled to gravity is described by an action of the form 

S = I d'x^ (^R - \d,^d>^v - \d,xd'x - ^(^, X)) , (29) 

where Mp is the reduced Planck mass, Mp = (SttG)^^/^. The homogeneous and isotropic 
FRW background with metric 

ds'^ = -dt^ + a(t)2rfx^ (30) 

is governed by the equations of motion (EOM) for the two scalar fields: 

(p + ^Hip + 1/^ = (31) 
X + 3i/x + Kx = 0. (32) 

Subscripts ip and x denote partial derivatives with respect to the corresponding field and a 
dot denotes a derivative with respect to the cosmic time, t. The gravitational background 
evolves according to Friedmann-Lemaitre equations: 

where H is the Hubble parameter, defined as = ^. 

To study the linear perturbations for this theory, we start with the longitudinal gauge for 
the metric [13] . In the absence of any anisotropic stress-energy tensor the scalar perturbation 
of the gravitational background reads: 

rfs2 = -(1 + 2<l>(t, x))rft2 + a{if (1 - 2<l>(t, x))dx2. (35) 

The scalar fields are also perturbed around their homogeneous parts, 

V9(t,x) = v2(t) + (5v9(t,x) and x(t, x) = x(t) + 5x(t, x), (36) 

These perturbations introduce an x dependence which was not present in the homogeneous 
and isotropic gravitational and scalar field backgrounds. To determine the perturbations, 
one therefore has to insert the perturbed metric and scalar fields into the full Einstein field 
equations and/or Bianchi identities, and the full scalar field EOMs. 

Since the perturbations of the metric and the scalar fields are not independent, it is useful 
to introduce gauge-invariant Mukhanov- Sasaki variables 

g^ = 5(^ + |:$ and Q^ = 5x+j^^. (37) 
18 



They represent the scalar field fiuctuations in the fiat gauge. It follows from the EOMs that 
their Fourier-components obey the coupled differential equationj^ 

+ 3HQ^ + + ^xx) Qx + C^^Q^ = , (39) 
with the following background-dependent coefficients 



_ 3ipx ^X" ^Vx I 



Ml 2M4ff2 2M^i72 + M^i/ + '^^^ ^^^^ 



Ml 2M4i72 2M^m^ MlH^ MlH^ ''''' ^ ' 

Following [H], we decompose the perturbations along and perpendicular to the trajectory in 
the (homogeneous) field space. The projection parallel to the trajectory is called the instan- 
taneous curvature, or adiabatic, perturbation whereas the one orthogonal to the trajectory 
is termed the instantaneous isocurvature, or entropy, perturbation. The velocity in the field 



space is 0" = + and we can define the polar angle in the field space as 

cos9 = Lp/a (44) 

It is now useful to define the following Mukhanov- Sasaki variable: 

Qo- = cos 61 Q,^ + sin 61 . (45) 

In the fiat gauge, Q^j represent the field perturbations along the velocity in the field space. Q^- 
is also related to the commonly used curvature perturbation, TZ, of the comoving hypersurface 
via 

n = -Q„ . (46) 
a 

Similarly the isocurvature perturbation is: 

5s = - sin 6Q^ + cos 6 . (47) 



^Even though, we will work with Fourier components hereafter, we will not show explicitly the subscript 
k which would denote the fluctuation with comoving wave-number k. 
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It describes field perturbation perpendicular to the field velocity in the field space and, by 
analogy with TZ, we can define a rescaled entropy perturbation, iS, through 

S = —5s . (48) 
a 

The transformations described above basically amount to introducing a new orthonormal 
basis in the field space, defined by vectors 

E^ = (E^,E^) = (cos0,sin0), (49) 
E, = {Ef, Ef) = (- sin e, cos 6) , (50) 

which turn out to be useful to express various derivatives of the potential with respect to 
the curvature and isocurvature perturbations. Employing an implicit summation over the 
indices J, J G {ip, x}, one thus finds 

K = EiVi , Vs = EiVi , (51) 

and 

= EiEiVij , = EiEiVij , Vss = EiEiVij . (52) 

for the first and second derivatives. 

By combining the Klein-Gordon equations for the background scalar fields, eqs. (l3Ti) and 
( 132!) . one obtains the background EOMs along the curvature and isocurvature directions 

— + ?,Ha + K = 0, (53) 

e = --. (54) 
(J 

With help of these one can show that the EOMs for curvature and isocurvature perturbations 
become 

+ ?>HQ, + [ +C„AQa + —Ss + C^sSs = 0, (55) 

. \ 91/ . 

5s + m5s +(— + Css]5s- + CsaQa = 0, (56) 



a? J a 
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with coefficients given by 

C -V -r^V + ^ + ^-^ f57) 

a. = 6i^- + + + ^ (58) 

= - ^ (59) 

A solution to these two coupled differential equations determines the metric perturbation $, 
which, in longitudinal gauge, is related to the comoving energy density 

+ \^3H + ^ j &Q„ + KQ<x + 2VJs (61) 

via the Poisson-like relation 

The power spectra of curvature (adiabatic) and isocurvature (entropy) perturbations are 
defined, respectively, as 

-PAk) = ^ (Q^kQ^k') S'{k - k') , Vssik) = ^ {6si6s^,) 6'{k - k') . (63) 

The curvature and isocurvature perturbations are then evolved by assuming initially, at 
conformal time Xj, a Bunch-Davies vacuum. Therefore, when the wavelength of the two 
types of perturbations is initially much smaller than the Hubble radius, k ^ aH, we impose 
the initial conditions 

^'^(^^^ = / X ATT ' = , , /TTT- • (64) 

a[Ti)y Ik a[Ti)y Ik 

Inside the horizon these two modes are independent, because their corresponding EOMs, 
eqs. (155|) and (156|1 . are independent in the limit k ^ aH. However, as we will see in detail 
later, this does not hold when the modes leave the horizon [l5l HI] . 
Finally, let us introduce the two field slow-roll parameters 

Vu = ^ (66) 

and 

e = e^^ + = . (67) 
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